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Based on the field theoretic formulation of the general epidemic process we study logarithmic correc- 
tions to scaling in dynamic isotropic percolation at the upper critical dimension d — 6. Employing 
renormalization group methods we determine these corrections for some of the most interesting time 
dependent observables in dynamic percolation at the critical point up to and including the next to 
leading correction. For clusters emanating from a local seed at the origin we calculate the number 
of active sites, the radius of gyration as well as the survival probability. 
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I. INTRODUCTION 

Spreading phenomena occur in nature in many kinds 
with examples ranging from epidemics or forest fires 
H EL 01 over the growth of populations H @ and 
the activity of catalyzers 0, H] to the formation of stars 
and galaxies 0- In general terms the spreading of non- 
conserved agents has the following scenario: An agent 
(e.g., an active site of a lattice, an infected individual 
or a burning tree) randomly activates one or more of its 
neighbors. In the next timestep these infected neighbors 
act as agents themselves and so on. Via this elemen- 
tary reaction, the activation spreads out diffusively in d- 
dimensional space. Competition between the agents for 
the resources of new activations limits their local density. 
Moreover, an agent become spontaneously deactivated 
after some time. The long term behavior of the process, 
assumed to be emanating from a seed or germ at the 
origin, depends crucially on the difference r between the 
deactivation and the activation rate. For r > r c the pro- 
cess spreads, approaching a homogeneous steady state, 
over the entire space. For r < r c the spreading will fi- 
nally ebb away as an inactive extinct state is approached 
(with possibly a static disturbance of the initial state, but 
only in a finite volume). The critical point r = t c sep- 
arates this endemic absorbing phase from the epidemic 
active phase. Spreading near r = r c constitutes a critical 
phenomenon and is described by universal scaling laws. 

There exist two fundamental universality classes of 
critical spreading phenomena depending on the nature 
of the debris of the elementary deactivation process. If 
the deactivated agents can recover so that they may be- 
come newly activated (so called simple epidemic), the 
process belongs to the directed percolation (DP) univer- 
sality class [l(| (for a review see, e.g., Ref. Here, in 
the active phase, one can have an epidemic surviving in 
loco. On the other hand, if the debris stays inactive for- 
ever (so called general epidemic), the spreading process 
becomes locally extinct. In this case, the process belongs 
to the dynamic isotropic percolation (dIP) universality 
class and the statistics of the clusters formed by the de- 



bris are described by the usual percolation theory (for 
reviews see, e.g., Refs. 0,0)- Here, the epidemic can- 
not, of course, survive in loco, but an infinite epidemic 
is nevertheless possible in the form of a solitary wave of 
activity. When starting from a punctual seed, this leads 
to annular growth (e.g., fairy rings in two dimensions). 

Most significantly perhaps, numerical simulations and 
renormalization group methods have contributed to our 
present understanding of spreading phenomena. Un- 
til recently accurate numerical investigations have been 
limited to lower spatial dimensions. However, due to 
the staggering pace of hardware improvements and the 
development of sophisticated algorithms, numerical re- 
sults for s p rea ding in high dimensions become available 
El 111 111 El Ell- Some of these results El El clearly 
indicate the importance of logarithmic corrections to 
scaling in the upper critical dimensions d = 4 for DP and 
d = 6 for dIP. Analytic work on logarithmic corrections to 
percolation started not long after renormalization group 
methods became available. Essam et al. E3 determined 
the logarithmic corrections for the probability of be- 
longing to an infinite cluster, the clusters' mean-square 
size S and the correlation length £ as functions of the 
deviation from criticality t — t c as well as the logarith- 
mic corrections at criticality to a "ghost" field H as a 
function of Poo. Aharony [2(j investigated logarithmic 
corrections in the context of universal amplitude ratios. 
More recently, Ruiz-Lorenzo [2lJ presented the logarith- 
mic corrections to S and £ at criticality as functions of 
the system size. In an upshot one can say that the previ- 
ous results on logarithmic corrections in percolation are 
(i) restricted to static percolation, and (ii) limited to the 
leading correction. 

For another system prominent in statistical physics, 
viz. linear polymers, logarithmic corrections observed 
in simulations have been very successfully described by 
renormalized field theory [22Ll23l |. It turned out, however, 
that the leading logartihmic corrections are not sufficient 
to yield a satisfactory agreement with the numerical data. 
To the contrary, it was found that it is crucial to include 
the next to leading correction. We expect a similar im- 
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portance of the second order correction in the percolation 
problem. 

The goal of this paper is to present analytical results 
for dIP in d = 6 that are, with reasonable expectation, 
accurate enough to to yield a satisfactory agreement with 
numerical simulations. We derive these results for some 
of the most interesting observables in dynamic percola- 
tion, namely the number N(t) of agents (active infected 
sites) at time t generated by a seed at the space- and 
time-wise origin (x = 0, t = 0), the survival probabil- 
ity P(t) of the corresponding cluster as well as the mean 
distance R(t) of the agents from the origin (radius of gy- 
ration) . 

The logarithmic corrections that we study here will 
not be found in real physical systems because the upper 
critical dimension 6 for dIP does not coincide with phys- 
ical dimensions, d = 2 or d = 3. However, our results 
are sure to be valuable with respect to numerical simu- 
lations. Our final analytic expressions are well suitable 
for comparison to numerical data. Moreover, our results 
define a non-universal time scale that signals the onset 
of asymptotic behavior. This time scale may be used 
to assess the effective significance of a given microscopic 
simulation model for the dIP universality class. 

Complementary to the work presented here we have 
investigated logarithmic corrections in DP. Our results 
on observables akin to the quantities studied here are 
certainly equally interesting and will be presented in the 
near future [24[ . Taking a third route, we explored loga- 
rithmic corrections in static IP with emphasis on trans- 
port properties. A paper on this subject |25| will be 
available soon. 

The outline of the present paper is the following: In 
Sec. ^ we briefly review the renormalized field theory of 
the general epidemic process (GEP) and previous renor- 
malization group results on dIP. Moreover, we conduct 
some general considerations about logarithmic correc- 
tions in the given context. Section II I II hosts the core 
of our analysis and contains the main results. Section llVl 
concludes the main part of this paper with a discussion 
of our results and several remarks. Details of our dia- 
grammatic perturbation calculation are relegated to an 
Appendix. 



II. RENORMALIZED FIELD THEORY OF THE 
GEP - A BRIEF REVIEW AND GENERAL 
CONSIDERATIONS ON LOGARITHMIC 
CORRECTIONS 

In this section we briefly review the field theoretic de- 
scription of the GEP and its renormalization. The aim is 
to provide the reader with background and to establish 
notation as well as known results that we need as we go 
along. Furthermore, we outline the general structure of 
the sought after logarithmic corrections. 

The GEP 0,1 i is a kinetic growth model that was 
introduced by Grassberger j2|J as a lattice model of dIP. 



On a coarse grained scale, the GEP can be minimally 
formulated by means of the Langevin equation in the Ito 
sense H3 

A _1 ?i(x, t) = V 2 n(x, t) — ™(x. t) — gn(x, t) 

x X f db'n(x,t') + C(x,t) , (2.1a) 

J — OO 

C(x,t)C(x',f) = A~V n(x, t)S(t - t')8{x - x') . (2.1b) 

Here n(x, t) is the density of infected particles at time 
t and space coordinate x. The variable r is essentially 
the rate difference mentioned in the introduction (shifted 
by r c ) and specifies the deviation from criticality. A rep- 
resents a kinetic coefficient. The Gaussian random field 
C(x, t) subsumes reaction noise and otherwise neglected 
microscopic details (the overbar indicates averaging over 
its distribution). Its correlation respects the existence of 
the absorbing state. Many other analytic terms are con- 
ceivable as contributing to Eqs. H2.ll) . These, however, 
turn out to be irrelevant in the sense of the renormaliza- 
tion group. Especially, a diffusional noise contribution, 
relevant for diffusion limited reactions with multiplicative 
noise, can be neglected. 

Langevin equations are only a convenient shorthand for 
a stochastic process. For the application of renormalized 
field theory, however, a path integral formulation of the 
GEP is more adequate than the Langevin equation (12.111 . 
The dynamic functional HHHHJ, or in a more recent 
termin olog y the response functional, of the GEP is given 
by IHUll 

J = J d d xdt\~s{\- 1 ^- t + {T-V 2 )+g{S-^S)^s. (2.2) 

In deriving J from the Langevin equation one exploits 
a rescaling form invariance of the response functional J 
that allows to equate g and g' . s(x, t) is proportional 
to n(x, t). s(x, t) is the response field corresponding to 
s(x, t). The S in the functional l|2.2|l stands for the den- 
sity of debris S(x, t) = A dt' s(x, t'). Note that 

s(x,t)< — >-S(x,-t) (2.3) 

is a symmetry transformation of the response func- 
tional |27| . 

J presents a vantage point for a systematic perturba- 
tion calculation in the coupling constant g. Most eco- 
nomically this calculation can be done by using dimen- 
sional regularization and minimal subtraction. Using this 
scheme, the critical point value r = r c is formally set 
to zero by the perturbational expansion. Generally, t c 
is a non-analytical function of the coupling constant g. 
Thus, we implicitly make the additive renormalization 
t — t c — > t. For background on these methods we refer 
to Refs. |33 . l33j |. An appropriate renormalization scheme 
is 

s->s = Z 1/2 s, 5^§=Z 1/2 §, (2.4a) 
A -> A = Z- 1/2 Z 1/2 X, t -> f = Z- x Z t t , (2.4b) 
g->g = Z-V 2 Z 1 J 2 g, G e g 2 ^u^. (2.4c) 
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Here, the ° indicates unrenormalized quantities, s — 6 — d 
measures the deviation from the upper critical dimen- 
sion. The factor G e = T(l + e/2)/(47r) d / 2 is introduced 
exclusively for later convenience, fj, is an external inverse 
length scale. Note that the renormalizations (|2.4I) pre- 
serve the invariance (|2.3|1 . The renormalization factors 
Z, Zr, and Z u are known to 3-loop order |3*4j . One of 
us [23 calculated Z to 2-loop order. In the following we 
will need the renormalization factors explicitly to 1-loop 
order to which they are given by 



Z = l-\ 
Z T = 1 
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(2.5a) 
(2.5b) 



where we abbreviated w = u(l). Solving this differential 
equation for e = 6 — d = yields 



I = l(w) = Ivw-^ltt 



exp 



1 

P 2 w 



(Pj - fob) 
Pi 



w+Q(w 2 



(2.10) 

where Iq is an integration constant. The remaining char- 
acteristics are all of the same structure, namely 



d]nQ(w) 
l^— = q{w) . 



(2.11) 



Here, Q is a placeholder for Z, Z, Z T , and Z\, respec- 
tively, q is an ambiguous letter for 7, 7, k, and respec- 
tively. Exploiting Id/dl = [3d/ dw we obtain the solution 



The critical behavior of the Green's functions G„,n = 
([s] n [s]™)( cm ") is governed by the Gell-Mann-Low renor- 
malization group equation (RGE) 

V^ + ^(n 1 + n^G n , h {{r,t};T,u;X,fi)=0 (2.6) 
with the differential operator 

-D„ = ^ + \(d x + TKd T + 0d a . (2.7) 

The Wilson functions appearing in the RGE are given to 
2-loop order bv pllpl 



4 /1895 ni o rl \u z 
7 = --.+ (^ r+ 91n3-51n4)- 



1 37 o 

7 = — u H u , 

' 6 216 ' 

5 193 2 

K = —U U , 

6 108 ' 



7-7 



7 /1747 

— u+ ( + 91n3-51n4 



12 V 54 



(2.8a) 
(2.8b) 
(2.8c) 



32 ' 
(2.8d) 



2" 72 V 2VS2 ' 93 ^ 3 - ) ) 16 

(2.1 



7 2 671 , / 414031 

p = — eu H — u u ' 



Q(w) = Qow 



qi/02 



exp 



/?2 2 



w + 0(w 2 ) 



(2.12) 

where Qo symbolizes a non universal integration con- 
stant. 

With the solutions to the characteristics the scaling 
behavior of the Green's functions is found to be 

G n> n({x.,t};r,w,X,n) 

= ( M o n(d+2)/2+ " ( ^ 2)/2 ^H n/2 ^H ?V2 (2.13) 

x G n ,n({lfiX, Z x {w)(l^) 2 Xt}; Z t {w)t/(hI) 2 ,w; 1, 1) . 

The flow parameter introduced via the characteristics is 
arbitrary. Thus, we have a freedom of choice that can be 
exploited to re-scale the relevant variables, viz. x, t and 
t -1 , so that they acquire a finite asymptotic value. For 
the goals pursuit in this paper, an appropriate choice is 



Z x (w)(lfi) 2 Xt = X . 



(2.14) 



where Xq is a constant of order unity. With this choice 
w and I tend to zero for \[i 2 t — > 00. Based on our 
choice (|2.14f) we introduce the convenient time variable 



7 . 



ln(t/*o) = -\n(t/t ) , 



(2.15) 



Note that we have stated to 3-loop order because the 
high order contribution improves our quantitative pre- 
dictions. In the remainder we will adopt a convenient 
abbreviated notation for the Wilson functions of the type 

f(u) = fa + fiu + f 2 u 2 H with / standing for 7, 7, 

K, C, and P, respectively. The meaning of the coefficients 
fo, fi and so on should be evident. 

The RGE can be solved by the method of characteris- 
tics. To this end one introduces a flow parameter I and 
sets up characteristic equations that describe how the 
scaling parameters transform if the external momentum 
scale is changed by varying p,(l) — The characteristic 
for the dimensionless coupling constant u reads 



dw 
~df 



P(w) 



(2.9) 



where to ^ X is a non universal time constant. From 
Eq. I|2.10|l and Eq. I|2.12|l . specialized to Z\, we get 



a\ In w + a 2 u> + 0(w 2 ) 



(2.16) 



for the derived time variable. The constants a\ and a 2 
are given by 



a 2 



P2C1 - 203 



1195 



2.37103, 



2/3 2 504 

Cl/?3-C2/3 2 P2P4-P! 

2/3 2 p\ 
1766273 101n2-91n3 279C(3) 



1016064 



64 



56 



(2.17a) 

(2.17b) 

= 7.68098 . 
(2.17c) 
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Using Eq. I|2.16[) we obtain for the dimensionless coupling 
constant as a function of time the expression 



s 1 exp 



hi.' 



fli- 



In s In s 1 



(2.18) 



Exploiting Eqs. ffTEty . (|2~TUj) . (|2~T2l) and l(2~T%)) we find 
that the observables to be considered are of the form 



A = Aq exp 



c In s + c' /In s Ins 1 
as+OlnsH hGq — 5—, —5-, -17 



(2.19) 

^4o is, like to, a non universal constant, a, fe, c and c' 
are universal numbers, a stems from mean field theory. 
b and c represent 1- and 2-loop renormalization group 
results, respectively, c' comprises contributions from the 
Wilson functions to 2-loop order as well as an amplitude 
to be determined in an explicit 1-loop calculation of A. 
This amplitude depends on Xq, as do s and Ao- Over 
all, a variation of X leaves c' invariant. 



III. LOGARITHMIC CORRECTIONS FOR THE 
OBSERVABLES OF INTEREST 

Equipped with important intermediate results as well 
as some knowledge of the structure of the Green's func- 
tions, we next determine the sought after logarithmic cor- 
rections. Since we already know the general form of the 
results, this part will be fairly brief. 



A. Number of active particles 

The number of active particles generated by a seed at 
the origin is given at criticality, t = 0, by 



N(t) = J d d xG h i(x,t;0,u; A,/i) 
= (Z(w)Z(w)) 1/2 [ d d x{iil) d 



x Gi,i(^r,Z A (w)(^) 2 At;0,w;l,l) (3.1) 
= (ZHZH) 1/2 G 1 , 1 (q = 0, X a ; 0, w; 1, 1) . 

Specializing solution l|2.12|l to Q = Z and Q = Z we 
obtain 



(71 \7( W 1 2 [(7i+7i) 
(Z(w)Z{w)) cxexp — — 

^P2 

' ( 72 + 12)132 ~ (71 + 71)03 

2/?I 



In w 



w + 0{w 2 ) 



(3.2) 



A perturbation expansion of the Green's function G\ t i 
brings about an amplitude Aj^(Xq) that we define via 

G M (q = 0, X ; 0, w; 1, 1) oc (l + A N (X )w + 0(w 2 )) . 

(3.3) 



This amplitude follows from our 1-loop calculation pre- 
sented in Appendix ITl as 



. . , 6 / 6 In 2 
A N (X ) = -(Z+-- — 



(3.4) 



where we used the shorthand notation 2 = lnA"o + C'e 
with Ce being Euler's constant. Collecting Eqs. (|3.1|l . 
(E3 and (E3J) we find 

N{t) = N (w- 1 + B N ) aN exp(c N w + 0{w 2 )) 



N'Js + B 



6 at In s + cn 



O 



In s In s 1 



(3.5) 



where N is a non universal constant, iVg is a non uni- 
versal constant slightly different from No, and Bn = 
Am /cln- The first row of 13.5fl and the result l|2.16|l con- 
stitute a parametric representation of the tuple (N, s) 
that is suitable for comparison to numerical simulations. 
The second row of (|3.5f) shows the more traditional form. 
The constants a at, &at, cat, and Bn are given by 



71 +7: 3 
a N = -l^ = -= 0.214286, 



(3.6a) 



on = ajv — = — norf> = —0.508078 , (3.6b) 



202 



2352 



72 +72 a 7l + 7l 

365 91n3-51n4 



= -0.206387, (3.6c) 



1568 112 

7, 3 ln2, 
B N = - {Z + - - — ) = 1-75Z + 2.22066 . (3.6d) 

Note from Eqs. (|2.15l) and (|3.5|1 that the arbitrary con- 
stant Z could be eliminated by a rescaling of the nonuni- 
versal time constant £o- This finding will also apply to 
the remaining results stated below. 

At this point we would like to warn against attempts 
to deduce the logarithmic corrections for dynamic quanti- 
ties from the logarithmic corrections calculated for static 
percolation. Grassber ger J15| . for example, exploited the 
results of Essam et al [T3 via replacing r by 1/i on the 
grounds that the critical exponent v t for the correlation 
time is one in mean-field theory. This reasoning leads to 
N(t) ~ [ln(t)] 2 / 7 [13. From Eqs. (^3 and how- 
ever, we see that the correct result, to leading order, is 
N(t) ~ [ln(t/to)] 3 / 14 ■ By merely using the mean field 
relation between r and t one misses contributions to the 
leading logarithmic term stemming from renormalization 
factors including Z [cf. Eq. I|3.2|l ]. Since Z is absent in 
static percolation |3a ] , one cannot deduce the logarithmic 
behavior of the dynamic quantity N(t) from the known 
results for static percolation. Likewise, it should not be 
attempted to combine our dynamic results, e.g. those 
for N(t) and R(t), to obtain predictions on logarithmic 
corrections in static percolation. 
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B. Radius of gyration 



by 



The mean square distance from the origin of the active 
particles is defined as 



R{tf 



J tfWG M (x,t) <9hiG M (q,i) 



2djd d xG ltl (x,t) dq 2 

(3.7) 

From the scaling form (|2.13|) it follows for r = that 



91nG 1 , 1 (q,i) 



dq 2 



S In G 1 ,i((//i)- 1 q,^xH(J/i) 2 At;0, u;;l,l) 



q=0 



= M 



dq 2 

_ 2 dlnGi ! i(q,X ;0,w;l, 1) 



<9<7 2 



q=0 

(3.8) 



q=0 



Incorporating the solutions to the appropriate character- 
istics and the results of Appendix HI 



a 



lnGi,i(q, \fj?t = X ; 0, w; 1, T 



with 



9g 2 

X (l + A fl (X ) w + O( w 2 )), 



, s 7 / 2 In 2 



q=0 



(3.9) 



(3.10) 



we find 



t" 1 ^ = RHw- 1 + B R ) aR exp(c R w + 0{w 2 )) 



bp In s + 



In s In s 1 



(3.11) 



with i?g and i?Q 2 being non-universal amplitudes. Here 
the constants ap, bp, cr, and Bp = Ap/ap are given by 



Ci 1 

ap = — — = - = 0.166666 . 
#> 6 



(3.12a) 



or = a-R 7, = = —0.395172 , (3.12b 

2/3 2 3024 ' v ' 

C2 Cl/?3 
CR ~ ~R W 

937 91n3-51n4 



= -0.128534, (3.12c) 
6048 112 ' v ' 

B R = 7 -{Z - H - = 1.752 - 1.33995 . (3.12d) 



C. Survival probability 

As shown in Ref. 37], the survival probability of an 
active cluster emanating from a seed at the origin is given 



P(t) = - lim (e- kAr S{-t)) 

k — ►QO 



(3.13) 



where Af = J d d x s(x, 0). For the purpose of actual cal- 
culations, it is convenient to rewrite Eq. I|3.13[) as 

P(t) = - lim (s(-t)) k 

k — >oo 

= — Go,i(— t, t, k — 00, u; A, fx) (3-14) 

where (•••)fe stands for averaging with respect to the 
response functional Ju that is obtained upon augmenting 
the original reponse functional l|2.2|l by a source k(t) = 
kd(t) conjugate to the field s: 



Jk = J+ / dtk(t)N(t) 



(3.15) 



With this source present, one no longer has (s) — 0. To 
avoid tadpoles in our perturbation calculation, we per- 
form a shift s — > s + M so that (s) — is restored This 
procedure leads to the new response functional 



Jk = J d d xdt (\s( 
+ (t - gM 
+ XgMsS -] 



dt 



1 



v 2 ) + g(s-^§))s 

-M+ XtM - y M 2 + k]s 



(3.16) 



Based on this functional we calculate Go,i = M to 1-loop 
order. Some details of this calculation are in Appendix^] 
We obtain 

G ,i(-X ,0, k = 00, to; 1, 1) 
oc vj- 1/2 {\ + A P (X )w + 0(w 2 )) . (3.17) 

with the amplitude Ap(Xq) reading 

A P (X )= 5 -(Z+1- 1 -^). (3.18) 

Recalling the scaling form (|2.13|) and our choice for the 
flow parameter we deduce that, for r = 0, 

P(t) = -Ziw^^lfGoA-Xo, 0, 00, w; 1, 1) . (3.19) 

Collecting, we then obtain 

tP(t) = P (ur 1 + B P ) ap exp(c P w + 0{w 2 )) 

bp In s + cp 



P^s + Bp) 



O 



'In s Ins 1 

. s 2 ' s 2 ' s 2 



(3.20) 



Po and Pq are simply related non universal amplitudes. 
The constants ap, bp, cp, and Bp = Ap/ap are given 
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by 



-P- Kl \ P r^ =Ta - 0-357143, (3.21a) 
2/j 2 14 

, 2/3 3 — 5975 „ 

bp = ap — = = -0.846797, (3.21b) 



2/3 2 



7056 



72 - 2C 2 . fl 2Ci - 71 
1637 91n3-51n4 



14112 



112 



= 0.0896074, (3.21c) 



Bp = 7 -(Z + 1- ^|^) = 1.752 - 0.918617. (3.21d) 



IV. DISCUSSION OF RESULTS AND 
CONCLUDING REMARKS 

As far as time dependent observables in percolation 
are concerned, we are not aware of any previous analytic 
work addressing logarithmic corrections. More general, 
we do not know of any work that has determined loga- 
rithmic corrections in percolation (static or dynamic, IP 
or DP) beyond the leading corrections. Here we went 
beyond the leading terms, and hence we are confident 
that our results compare well with simulations, perhaps 
even quantitatively. For linear polymers it turned out 
that the knowledge of the leading logarithmic correction 
is not sufficient for a good agreement between simulation 
data and theory. Rather, the next to leading corrections 
turned out to be crucial in comparing numerical and an- 
alytical results. We expect the same for percolation. In- 
deed, preliminary Monte Carlo results corroborate this 
expectation 0,0. 

Our results define a non-universal time-scale to. For 
times t greater then to we expect the validity of our 
asymptotic expansions. The time-scale to can be utilized 
as a measure of quality for different microscopic models 
of dynamical percolation. Thus, our results may guide 
those performing simulations in choosing the most effi- 
cient model. 

It is interesting to note that the time-scale to has 
an analog in quantum chromodynamics. For times 
greater than to the model becomes asymptotically free. 
Thus, with the exchange of an infrared(IR)-free the- 
ory to an ultraviolet (UV)-free theory, to corresponds 
to the hadronization scale of quantum chromodynam- 
ics. The dependence of our results on this dimensional 
non-universal parameter to parallels therefore the phe- 
nomenon of dimensional transmutation in renormaliz- 
able asymptotically free quantum field theories that are 
naively scale-free. 

Our results feature a mutual non- universal constant, 
viz. Z. This constant could be eliminated by rescaling 
to- One might be tempted to think that one could elim- 
inate the entire amplitudes, An(Xq) and so on, from 
our results by rescaling to, and that the calculation of 
the amplitudes is hence superfluous. One has to keep 
in mind however, that one has to choose to consistently 



for all observables. Thus, one cannot remove the ampli- 
tudes simultaneously from all the observables, and their 
calculation is indeed necessary. 

We refrain from eliminating Z from our results because 
it might be exploited, due to its non-universality, as a fit 
parameter. By fitting Z one can compensate partially for 
the effect of higher order terms that have been neglected 
in our calculations. In this sense one can think of Z as 
mimicking these higher order terms. 

When written as an explicit function of time the ob- 
servables of interest have fairly complicated formulae. 
Using the parametric representation in terms of the ef- 
fective coupling constant w eases this situation. More- 
over, the time and the observables posses a nicely sys- 
tematic expansion in w so that it is straightforward to 
keep track of the different orders in perturbation the- 
ory. In the traditional from, the orders are not so clear 
cut because nested functions of logarithms have to be 
compared. The parametric representations can be con- 
veniently compared to simulations. Essentially, one just 
needs to make parametric plots of the tuples (TV, s), (R, s) 
and (P, s) and then compare the numerical data to these 
plots. 

In order to improve the accuracy of our results, one 
needs a refined quantitative knowledge of the Wilson 
functions. Whereas the e expansion results for the criti- 
cal percolation exponents have been improved by resum- 
mation techniques such as Pade-Borel resummation [34| , 
this kind of refinement has not yet been achieved for the 
percolation Wilson functions. Here lies an opportunity 
for useful future work for the ambitious. Another possi- 
bility for future work is to improve the results on static 
percolation mentioned in the introduction by calculating 
the next to leading logarithmic corrections. With the 
kind of field theoretic methods that we applied here, this 
is a reasonable task. 

Apparently, firm numerical results that are suitable for 
comparison to our analytical results are not available yet. 
We hope, however, that our work triggers increasing ef- 
forts in this direction, and that corresponding numerical 
results will be available in the near future. 
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APPENDIX: EXPLICIT CALCULATION OF 
GREEN'S FUNCTIONS 

In this Appendix we outline our 1-loop calculation of 
scaling functions belonging to the Green's functions G\ t \ 
and Go,i. In particular, we compute the amplitudes 



FIG. 1: Self-energy E(q,t) at 1-loop order. 



An{Xq), Ap(X ), and Ar(X ) entering the logarithmic 
corrections. 



1. The Green's function Gi,i 

A hrst step of any diagrammatic perturbation calcula- 
tion is, of course, the determination of the constituting 
elements. From the response functional (|2.2|) we gather 
the Gaussian propagator 



G(q,t) =6(t) exp[- X(r + q 2 )t] 



(A.l) 



and the 3-leg vertices Xg and —X 2 g9(t — t'), where 9(t) 
denotes the step function. With these elements, the self- 
energy S(q, t) is given at 1-loop order by the diagram 
depicted in Fig. ^ This diagram stands for the mathe- 
matical formula, 

E(q, t) = -X 3 g 2 f dt' f G(p, t')G(q -p,t). (A.2) 
Jo Jp 

After integrating out the loop momentum we can rewrite 
E(q, t) as 



E(q,i) 



(A 5 ) 2 



■(Xt) 



l-d/2 



ds 



(4 7r )d/2 

exp[-At((l + s)t + sg 2 /(l + s))} 
(1 + s) d ' 2 



(A.3) 



For our purposes we need Green's or connected corre- 
lation functions rather than vertex functions. Hence, we 
have to consider Feynman diagrams with external legs 
attached rather than amputated diagrams. The Green's 
function G% t x is determined by the Dyson equation 

Gi,i(q,t) = G(q,t) + / dt'f dt" G(q,i - t') 
Jo Jo 

xE(q,f'-(")G(q,f") + -. (A.4) 

A diagrammatic representation of the Dyson equation is 
given in Fig. |21 Upon substituting l|A.3|) into (|A.4I) we 
obtain after an integration 



G 1 , 1 (q,t) = G(q,t) 



1 - 



i{Xn 2 t) 



2 f) e/2 



ds 



r(l + e/2)7 (! + «)«*/ 



x / dx(l-x)x 1 d / 2 exp(a(s)x) 
Jo 



(A.5) 



FIG. 2: Dyson equation (IA.4I to 1-loop order 

Here, we used the shorthand notation 

,2 



a( a )=( T i_-(l + a )r)At. 



(A.6) 



Now, we set r = and expand Eq. i|A.5|) to order q 2 . 
The integrations are easily performed. After e-expansion 
we get 



Gi,i(q,*)=G(q > *) 



u{X^i 2 t) 



2 f y/2 



r(l + e/2)VV4e 16 



(t 



9 In 2 



(. 



7 



7 In 2 



V12e 36 



24 



Xq 2 t 



(A.7) 



The next step is to remove the e poles by employing 
the renormalization scheme (|2.4I) . Letting Gi,i — > G14, 
A — > A, and using 

G M = (^)-V2 Gl)1 = (l - g)G M (A.8) 

as well as 

7« 



(A.9) 



we observe that the e poles are indeed removed. For the 
renormalized Green's function we obtain 

Gi,i(q, t) = 1 + ^ ( ln(A/A) + G £ + | - ^ 



g(ln(Vi)+C*-§-^) (A.10) 



Two results important for the logarithmic correction 
can be extracted from l|A.10|) . Upon setting q = we 
find 



G M (q = 0, XfiH = X q ;t = 0, w; 1, 1) 

3/ 3 In 2 
^ 1+ 8( Z+ 2-^'^ 



(A.ll) 



and hence the amplitude An(Xo) as stated in Eq. (|3.3|l . 

Moreover, we get 



- X- 1 A hiGi,i(q, X^t = Xq; t = 0, w; 1, 1) 

7 / 2 In2\ 
^ 1+ 2i( Z -3--) W ' 



9 2 =o 
(A.12) 



which leads to our result for A/{(Ao) given in Eq. I|3.10|l . 
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a) 



differential equation (|A.13(1 we find the modified Gaus- 
sian propagator 



FIG. 3: (a) The new vertex -X 2 gM 6(t-t') and (b) the 1-loop 
tadpoie diagram T(t). 



2. The Green's function Go,i 

Now we determine Gq i as required in (|3.14|) . The dia- 
grammatic elements associated with the functional 
comprise the two vertices encountered in Appendix^ In 
addition, there is a third vertex, viz. the one depicted in 
Fig. EH. The Gaussian propagator for the new functional 
has to be determined from the differential equation 



X^dt+r-gM^+q 2 G(<i,t,t') = \- 1 5(t-t'). (A.13) 



To avoid tadpoles, M(t) has to satisfy the differential 
equation 

M(t) - XrM(t) + y M(i) 2 - k{t) + T(t) = . (A. 14) 

At 1-loop order, the tadpole T(t) is given by the diagram 
shown in Fig. EJd. 

The initial and terminal conditions for the fields neces- 
sitate the ansatz M(t) = -6>(-i)if (-i)" 1 . The type of 
the source term, k(t) = kS(t) with k — > oo, demands the 
initial condition K(0) = 0. With this information, the 
differential equation (|A.14|) can be transformed without 
much effort into the integral equation 



K(t) 



,Art 



2r 



dt'e- XTt 'K{t') 2 T{-t') + ^j 



(A.15) 

At mean field level, the solution to Eq. (|A.15|I is given 
by 



K o(t) = ^(e XTt -l). (A.16) 
Inserting the corresponding Mo(t) = —Ko(-t)^ 1 into the 



(A- 17) 

Having the modified Gaussian propagator at our de- 
mand, we are now in the position to calculate the dia- 
gram depicted in Fig. EJa. Eventually, we obtain 

K(t) 2 T(-t) = X 3 g 2 K (t)- 2 f dt' f dt" 
K {t')K n {t") 2 exp[Ar(2i -if - if')] 



[4ttA(2< -t' - t")] d / 2 



(A.18) 



The further evaluation of Eq. (|A.18|) is fairly straightfor- 
ward for t = 0. Away from criticality, the calculation is 
more challenging and will be addressed in a future publi- 
cation |3J|. Here, we find for r = and after e expansion 



, . Xg 3 (Xt) £ / 2 ( 5 5 llln2\ , , , 

K(t) 2 T{-t) = - y fj }> - + - - — — . (A.19) 



(47r) d / 2 V8e 8 16 
Insertion of this intermediate result into l|A.15(l yields 



Kit) = & 



. u{Xn 2 t) £ l 2 / 5 + 5 11 In 2 



r(l + e/2) \Ae 



(A.20) 



Next, we renormalize. Indicating the consistency of our 
previous steps, the appropriate combination of renormal- 
ization factors (ZZ jZu)^ 1 ! 2 = 1 + 5u/(4e) H — • cancels 
the £ pole in (|A.20I) . The renormalized K(t) reads 



K(t) 



gXt 



5u 



(]n(\(J?t)+C E + l 



llln2 



8 V v r ' " 5 

(A.21) 

Exploiting Go,x(—t) = K(t) 1 and X/j, 2 t = Xq as well as 
recalling the definition of Z we finally obtain 



(4tt) 3 / 2 X 



2 

,,-1/2 



G ,i(-A/i 2 t = X ;t = 0,w;l,l) 



, 5/ llln2 
1+- Z+l — )w 



(A.22) 
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